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ABSTRACT 

The contribution of weakly-magnetized (B ~ 10 9 G) neutron stars accreting the 
interstellar medium to the diffuse X-ray emission observed in the Galactic Center is 
investigated. It is shown that, under rather conservative assumptions about the neutron 
stars and gas distributions, the accretion luminosity can account for a sizable fraction, 
possibly most, of the detected X-ray flux in the 2.5-7 keV band. In particular, model 
results are compared with Granat data and show a general agreement in both the flux 
energy and radial distributions. 



Subject headings: Accretion, accretion disks — Galaxy: center — stars: neutron - 
X-rays: general — X-rays: stars 
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1. Introduction 

It was suggested long ago in a seminal paper by 
Ostriker, Rees, & Silk (1970) that old isolated neu- 
tron stars (ONSs) can be roused from their lethargy 
by the accretion of the interstellar medium, revealing 
themselves as weak, soft X-ray sources. About 20 
years later it was realized that the capabilities of the 
ROSAT satellite were unprecedented in matching the 
challenge of observing ONSs (Treves & Colpi 1991; 
Blaes & Madau 1993; Madau & Blaes 1994; Zanc et 
al. 1996). The basic criteria for discriminating ONSs 
from other classes of X-ray sources should be their 
very low luminosity, extreme X-ray to optical flux 
ratio and spectral softness. In accordance with the 
pioneering suggestion by Shvartsman (1971), giants 
molecular clouds were recognized as very favourable 
sites for the detection of accreting ONSs (Blaes & 
Madau 1993; Colpi, Campana, & Treves 1993). 

Despite the increasing efforts, up to now only two 
X-ray sources have been proposed as ONSs accret- 
ing the interstellar medium: MS 0317.7-6647 (Stocke 
et al. 1995) and RXJ185635-3754 (Walter, Wolk, & 
Neuhauser 1996). In both cases the arguments in 
favour of their identification are strong but not com- 
pletely compelling. Systematic searches on available 
X-ray data are under way (see e.g. Danner, Kulka- 
rni, & Hasinger 1994; Bclloni, Zampicri, & Cam- 
pana 1996), but the positive detection of ONSs, which 
could be as many as 10 9 in the Galaxy, appears still 
a difficult goal to achieve. 

Because of the intrinsic difficulties in detecting sin- 
gle objects and of the expected relative abundance of 
ONSs in the Galaxy, it is natural to consider their 
overall emission and more specifically their contribu- 
tion to the soft X-ray background, as already sug- 
gested by Ostiker, Rees, & Silk (1970). Assuming 
B ~ 10 9 G and N ONS ~ 10 9 , (Zane et al. 1995) have 
recently shown that ONSs accreting the interstellar 
gas can account for a non-negligible fraction of the 
observed excess in the soft X-ray background which 
is not extragalactic in origin but could be produced 
by a Galactic population (Maoz & Grindlay 1995). 

Since the diffuse emissivity of accreting ONSs de- 
pends on their spatial concentration and on the den- 
sity of the interstellar medium, it is of particular in- 
terest to consider the Galactic Center, a site where 
both the stellar and gas densities exceed by orders of 
magnitude those of other regions of the Galaxy. At 
the same time the Galactic Center is a well known 



source of diffuse X-ray emission, first detected by 
Uhuru (Kellogg et al. 1971) and then studied by vir- 
tually all X-ray missions. The possible association 
of the Galactic Center diffuse X-ray source with ac- 
creting ONSs was originally suggested by Maraschi, 
Treves, & Tarenghi (1973). 

In this paper first we review the X-ray observations 
of the Galactic Center. We then discuss the interstel- 
lar medium, the stellar and the expected neutron star 
distributions. The calculation of the emission due to 
accreting ONSs is presented and our results are com- 
pared with observational data. Discussion follows. 

2. The Galactic Center 
2.1. The X-ray Emission 

The Galactic Center is one of the more widely ex- 
plored regions of the sky in the X-rays and observa- 
tional efforts appear indeed motivated in the light of 
the complexity of the source. The inner 100 pc of the 
Galaxy exhibit, in fact, a region of diffuse emission 
together with a number of point-like sources. The 
strong absorption in the direction of the Galactic Cen- 
ter makes its appearance substantially dependent on 
the energy band (soft or medium X-rays) in which the 
source is observed. Here we briefly outline the current 
status of X-ray observations of the Galactic Center, 
focussing our attention on the diffuse component. 

The presence of a weak, diffuse emission from 
about ~ 1 deg 2 was already suspected in the Uhuru 
2-10 keV data (Kellogg et al. 1971) and first con- 
firmed by Einstein in the 0.5-4 keV band (Watson 
et al. 1981). Observations with 10' resolution in the 
2-15 keV range, performed with Spartan 1 (Kawai et 
al. 1988) and Spacelab 2 (Skinner et al. 1987; Skin- 
ner 1989), confirmed the diffuse emission at higher 
energies. Spectral measurements by Ginga, with an 
angular resolution of ~ 1°, have shown the presence 
of a strong emission line at 6.7 keV, which was iden- 
tified with the Ka line of He-like iron (Koyama et 
al. 1989; Yamauchi et al. 1990). The continuum 
appears rather flat and is well fitted by a thermal 
bremsstrahlung. A temperature of 12.7 ±0.4 keV was 
calculated on the basis of TTM observations by Not- 
tingham et al. (1993). Such temperature, however, 
is more than one order of magnitude too high for a 
plasma to be confined in the Galactic Centre potential 
well. 

The ART-P telescope on board Granat observed 
the Galactic Center region with 5' resolution in the 
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2.5-30 keV band with rather long exposure times 
(Sunyaev, Markevitch, & Pavlinsky 1993; Marke- 
vitch, Sunyaev, & Pavlinsky 1993). The intensity 
profiles obtained in 4 energy bands, after subtracting 
point-like sources, show the presence of an elliptical, 
extended source, as first suggested by Kawai et al. 
Spectra from the total 1° x 1.5° ellipse and from a 
central region 30' wide were produced, confirming a 
flat bremsstrahlung like continuum; a strong absorp- 
tion feature at 8-11 keV was also detected. Sunyaev, 
Markevitch, & Pavlinsky (1993) noted that the struc- 
ture of the diffuse emission differs substantially below 
and above ~ 8 keV. While the lower energy compo- 
nent is thermal and roughly elliptical, the hard emis- 
sion comes from an elongated region, parallel to the 
Galactic plane, that resembles the distribution of Gi- 
ant Molecular Clouds. This led Sunyaev, Markevitch, 
& Pavlinsky to the conclusion that the diffuse emis- 
sion consists of two components and that the high en- 
ergy portion of the spectrum may be due to Thomson 
scattering of hard photons on the dense material of 
the clouds. In this picture the bremsstrahlung tem- 
perature could be lower than the previous estimate 
(Markevitch, Sunyaev, & Pavlinsky give Tb rems ~ 3 
keV) , easing the problem of confining the hot gas in 
the Galactic Center. 

The best resolution X-ray map of the Galactic 
Center was obtained with ROSAT in the 0.8-2 keV 
range (Predehl & Triimpcr 1994). In order to explain 
the lack of X-ray sources at the position of Sgr A* , an 
interstellar absorption higher than 2 x 10 23 cm -2 was 
invoked by Predehl & Triimpcr. Preliminary reports 
of ASCA observations (Koyama et al. 1996) indicate 
the presence of several metal lines besides iron. In 
particular, a 6.4 keV fluorescent Ka component ap- 
pears superposed to the 6.7 keV emission feature. 

2.2. The Gas Distribution 

The structure of the interstellar medium (ISM) in 
the central 10 2 parsecs of the Galaxy has been ex- 
tensively investigated and a detailed review can be 
found in Genzcl, Hollenbach, & Townes (1994). In- 
formations about the dust distribution are obtained 
from the re-radiation of UV and visible photons into 
the infrared continuum, atomic and ionized compo- 
nents are observed directly in the 21 -cm line while 
the more abundant molecular gas is sampled by the 
millimetre, submillimctre and infrared lines of trace 
molecules (CO, OH, HCN, ...). The resulting picture 
shows that the Galactic Center is a region charac- 



terized by a strong concentration of dense interstel- 
lar material, with an average density of gas and dust 
10 2 — 10 5 times higher than in the rest of the Galaxy. 
The central 5 parsecs contain a circum-nuclear disk 
(CND) of orbiting filaments and streamers with its 
inner edge at ~ 1.5 pc from the center. This disk is 
probably fed by the infalling gas from denser molec- 
ular clouds at r > 10 pc and drops streamers in the 
central region. The inner region (r <J 1.5 parsecs: 
the central cavity and the mini-spiral) are compara- 
tively devoid of material and the average gas density 
is, at least, one order af magnitude lower than in the 
CND. On larger scales, surveys of 2.6 mm CO and 
far-infrared dust emission with IRAS (Dame et al. 
1987, Deul & Burton 1988) show that - 10 s M Q of 
gas (~ 10% of the total Galactic ISM) are contained in 
the inner few hundred parsecs. Using a new CO-to— 
H 2 conversion factor, Sofue (1995a) estimated a value 
for the gas mass of - 4.6 x 10 7 Mq within w ~ 150 
pc, where w is the distance from the center in the 
plane of the sky. The corresponding average gas den- 
sity turns out to be ~ 100 cm , in agreement with 
the value reported by Genzel, Hollenbach, & Townes 
(1994) for the inner 100 pc. The total mass of gas and 
dust is only 1-10% of the stellar mass and does not 
contribute significantly to the Galactic gravitational 
potential. The dust comprises only about 1% of the 
ISM mass which is mainly in the form of molecular, 
atomic and ionized gas. A fraction of the interstellar 
material is organized in dense molecular clouds and 
macro-structures: within ~ 500 pc the filling factor 
is — 10% for clouds with n ~ 10 3 cm~ 3 and ~ 0.05% 
for clouds with n ~ 10 5 cm~ 3 (see, e.g., Campana & 
Mcrcghctti 1993 for a discussion). As recently dis- 
cussed by Sofue (1995a), the gas distribution in the 
central ~ 1 kpc region is dominated by a rotating ring 
with n w 10 3 cm~ 3 located at ~ 120 pc while the gas 
density drops by one order of magnitude outside. An 
expanding, roughly spheroidal, molecular shell is also 
observed at r ~ 180 pc, \z\ <J 50 pc, with the gas 
mainly concentrated at intermediate latitudes (Sofue 
1995b). 

Using IRAS observations of the 3° x 2° region 
around Sgr A, Cox & Laureijs (1989) proposed that 
the dependence of the gas mass on the projected dis- 
tance from the Galactic Center is 

m gas {< vj) ~ 2 x lO 3 ro 1 - 8 M . (1) 

Such distribution holds approximately up to 300 pc 
from the center and implies, on such scales, a nearly 
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constant surface density. As discussed by Cox & Lau- 
reijs, this result is very sensitive to the assumed dust 
temperature and the value of the total gas mass is 
correct to about 30 % within 300 pc (50 % at larger 
distances) . Clearly, the averaged volume density pro- 
file can be derived only by de-projection and any dis- 
tribution matching the constraint of constant surface 
density is a priori acceptable. In particular, if the ex- 
tension of the gas along the line of sight is roughly 
constant when w varies, it is not unreasonable to as- 
sume an homogeneous gas distribution with an aver- 
aged, constant value for n. Within 225 pc, the av- 
eraged value of n derived from IRAS data is ~ 60 
cm~ 3 , assuming an elliptical gas distribution with an 
axial ratio 0.7 and using a value of 3.6 x 1O 7 M0 for 
the total mass (Cox & Laureijs). However, in the in- 
ner ~ 100 pc this value is probably an underestimate 
(see the preceeding discussion) and n = 10 2 cm~ 3 
seems to be closer to observations. In the following 
sections we will focus our attention on the contribu- 
tion of accreting ONSs to the diffuse emission from 
the Galactic Center. The accretion rate will be es- 
timated assuming an homogenous distribution of the 
ISM with 

n= 100 cm" 3 nj<100pc (2) 

Outside this region, in order to correct the X-ray 
emission for interstellar absorption, we need an esti- 
mate of the gas density in the Galaxy. This is known 
to be of order unit and is a parameter in our model, 
adjusted in such a way to give a total column den- 
sity, integrated over 8.5 kpc at z = 0, always greater 
that N H = 2 x 10 23 cm~ 2 . As discussed by Predchl 
& Triimper (1994), such a large value of Nh appears 
compelling if the deficit of X-rays sources observed 
by ROSAT toward Sgr A* has to be accounted for. 

2.3. The Star Distribution 

Both observations and theory (see e.g. Bailey 1980; 
Allen, Hyland, & Hillier 1983; Sanders 1989; Morris 
1993) provide convincing evidence that in a very large 
range of Galactocentric distances, from few arcsec- 
onds to few degrees from the center, the star volume 
density scales approximately as r~ 2 (r is the galacto- 
centric spherical radius) in agreement with the pre- 
diction of the isothermal cluster model for the region 
outside the core radius a. The value of the core radius 
derived from observations ranges between ~ 0.04 and 
~ 0.8 pc (see e.g. Morris 1993; Genzel, Hollenbach, 
& Townes 1994 and references therein). Within the 



core radius, data derived from the dynamics of both 
the stellar and gas components agree quite well in sug- 
gesting that an additional dark mass of ~ 3 x 10 6 Mq 
is present in the central region. This mass could be 
in the form of either a single, massive black hole or, 
if the central density is not in excess of 10 8 Mq pc~ 3 , 
a cluster of many compact remnants (Morris 1993; 
Genzel, Hollenbach, & Townes 1994). According to 
this picture, the mass volume density of the stellar 
component is given by 

p(r)=p c r- L8 M pc- 3 , (3) 

where p c is the central density which can be inferred 
from the estimated mass inside the core radius; typi- 
cal values for p c are a few 10 5 Mq pc -3 (Sanders 1989; 
Morris 1993). In the following we assume that this 
profile adequately describes the stellar distribution in 
the region beyond 1 pc from the Galactic Center and 
we take p c = 4 x 10 5 Mq pc -3 . 

The fraction / of the total mass comprised in neu- 
tron star is a free parameter of our model and we 
assume that all neutron stars have the same mass, 
M ns = 1.4Mq. For a Salpeter initial mass function 
(IMF) 

N(M) oc M- a (4) 

with a = 2.35 and assuming that all stars with ini- 
tial masses between ~ IOMq and ~ 50Mq have left 
neutron stars, / turns out to be ~ l %. With a fixed 
a, the previous formula can be assumed to describe 
the IMF of an already formed galaxy for masses in 
a given range, so possible differences due to the evo- 
lution of the bulge in the early Galactic history are 
ignored. However, as suggested by Morris (1993), the 
initial mass function for stars formed in the Galactic 
Center may be quite different from that averaged over 
the disk of the Galaxy, due to the different physical 
conditions in the star-forming clouds. In particular, 
throughout most of the central region, strong tidal 
forces increase the value of the limiting density for 
a cloud to become self-gravitating and this, in turn, 
acts toward inhibiting star formation. This implies 
also that close to the center self-gravitating clouds 
are more and more unlikely to be present and that 
those in which the density is above the tidal limit are 
presumably supported by the strong magnetic pres- 
sure. Large internal velocity dispersions and higher 
temperatures within the clouds tend to increase the 
Jeans mass up to an order of magnitude with respect 
to the typical value in the rest of the Galaxy. In the 
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inner few parsecs the most effective way of forming 
stars seems to be external cloud compression due to 
collisions with other clouds and/or with shock fronts 
produced by supernova explosions and other forms 
of nuclear activity. During the compression the gas 
temperature and the value of the Jeans mass are in- 
creased. Basing on these arguments Loose, Kriigel, 
& Tutukov (1982) argued that the mass spectrum in 
the Galactic Center is probably skewed toward higher 
masses, with a lower cutoff at M m i n ~ 1 — 3Mq 
(M min = 0.08Mq is often assumed in the disk). In 
addition, the value of a itself turns out to be smaller 
than in the disk, so / = 0.01 should be regarded as a 
lower limit. Corrections due to the effects of dynami- 
cal friction over the lifetime of the Galaxy have been 
considered by Morris (1993), but they have not been 
included here since they are more relevant for black 
holes due to their larger masses. Black hole remnants 
are expected to dominate the central cluster and their 
total mass is a large fraction of the dynamically in- 
ferred one. On the other hand, since mass segregation 
is not expected to modify strongly the distribution of 
neutron stars remnants, in the following we will use 
a number density of neutron stars given, outside the 
core radius, by 



Since the total luminosity emitted in the accre- 
tion process depends on both the gas density and the 
star velocity, the velocity distribution f(v) of neu- 
tron stars is needed in order to estimate the cumu- 
lative emission associated with the NS population. 
The velocity distribution of Galactic ONSs can be 
obtained assuming a suitable distribution of NSs at 
birth in phase-space and following the evolution of a 
large number of orbits in the Galactic gravitational 
potential (see e.g. Blaes & Madau 1993, Zanc et al. 
1995 and references therein). This kind of approach 
proved very useful in calculating f(v) on large scales 
but, due to the particular conditions of the Galactic 
Center, results obtained from a detailed analysis are 
not necessarily more correct than a simpler estimate. 
As discussed by Sellgren et al. (1990, see also Genzel, 
Hollenbach, & Townes 1994), the velocity dispersion 
a v for low-mass stars in the central region is nearly 
constant between 0.6 pc and 100 pc, with a v ~ 75 
km/s. At very small r observational data indicate 
that a v increases up to ~ 125 km/s. Here we assume 
that NS remnants were able to reach energy equipar- 
tition with the field stars, establishing a Maxwellian 



distribution 

/(«) = ^cxp(-3x 2 /2) , (6) 

where x = v/vq, vq = a v and typical values for a v 
are in the range 75 — 125 km/s. For a typical mass 
of ~ IMq, it can be easily shown that the relaxation 
time t re iax (see e.g. Binney & Tremaine 1987) is 
<J 10 10 yr in the inner ~ 10 pc. Since t re iax is a 
measure of the time required for deviations from a 
Maxwellian distribution to be significantly decreased, 
in the innermost region our hypotesis is justified. En- 
ergy equipartition between populations of different 
mass is reached in about the same time. In addi- 
tion, the evaporation time turns out to be ;> 10 10 yr 
for r <; 10 pc, so the irreversible leakage of stars from 
the system due to stellar encounters is negligible over 
the Galaxy lifetime. A similar approach was used by 
Morris (1993) for the velocity distribution of the cen- 
tral cluster of black hole remnants. Extending the 
Maxwellian assumption up to ~ 100 pc is motivated 
on the basis of observational data. The total mass 
scales approximately as r here and, as previously dis- 
cussed, this is just what is expected in an isothermal 
cluster model. Actually, a lowered Maxwellian distri- 
bution is probably more reasonable, but the correc- 
tion in the high energy tail (above ~ v esc ~ 2a v ) is 
not very important as far as the fraction af accreting 
ONSs is concerned, and will be not included in our 
calculations. 

3. Results 

In this section we estimate the contribution of the 
neutron star population to the diffuse emission near 
the Galactic Center, using the NS and gas distribu- 
tions considered in the previous sections. For a gen- 
eral discussion of neutron stars accreting the ISM we 
refer to Treves & Colpi (1991) and Blaes & Madau 
(1993). 

The spectral properties of the radiation emitted 
by accreting ONSs have been recently investigated by 
Zampicri et al. (1995). They have shown that, even 
if the atmosphere is in LTE, the emergent spectrum 
is harder than a blackbody at the neutron star effec- 
tive temperature and the hardening ratio, typically 
<~ 1.5 — 3, increases with decreasing luminosity. This 
spectral hardening, produced by the decoupling of dif- 
ferent frequencies at different depths, is present inde- 
pendently of the assumed geometry of the accretion 
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flow. If ONSs retain a relic magnetic field which chan- 
nels accretion onto the polar caps, the spectrum is 
harder because of the reduced emitting area (see e.g. 
Treves & Colpi 1991) and the two effects add together, 
making the emission of medium X-rays (^4 — 5 keV) 
possible in the dense region of the Galactic Center. 
In the following we assume that the emitted spec- 
trum is that calculated by Zampicri et al. and that 
ONSs have a relic magnetic field B = 10 9 - 10 10 G. 
Note that in this case the luminosity released in the 
cyclotron line is negligible (Nelson et al. 1995). 

The monochromatic flux emitted by ONSs in a re- 
gion of volume V centered at r — can be calculated 
from the expression 

F v = / n (r) wdwdzdcf) x 
Jv 

J f^4^W dV er g cm ~ 2frl > ( 7 ) 

where L u is the monochromatic luminosity at the 
source corrected for the interstellar absorption. In the 
previous expression the volume element is expressed 
in terms of the cylindrical coordinates (zu, z,<p): z is 
along the line of sight and w has been defined in sec- 
tion 2. R = 8.5 kpc and R = (i?g+tu 2 -ZRuz+z 2 ) 1 / 2 
are the distance of the Galactic Center and of the star 
from the Sun. The integral in equation (7) has been 
evaluated for \z\ < z max and for different values of zu. 
IRAS data discussed in section 2 show that the dense 
material extends for a few hundred parsecs around 
the center, so we assume z max — 300 pc. The un- 
absorbed monochromatic flux at the source depends 
on the accretion rate which, in turn, is a function 
of v and n. To avoid the direct calculation of the 
spectrum for each value of the parameters within the 
multiple integral, we have computed a set of models 
for different values of M. The spectral distribution 
corresponding to any given pair of values of v and n 
is then obtained by spline interpolation. The total 
emitted flux and the total flux for square degree in a 
fixed energy band are calculated integrating equation 
(7) over frequencies. All multiple integrals have been 
evaluated numerically and the cross sections for the 
interstellar absorption has been taken from Morrison 
& McCammon (1983). 

The region within ~ 1 pc from the center (central 
cavity) , where our assumed star and gas distributions 
cease to be valid, has been excluded from the integra- 
tion domain in all models. The X-ray emission from 
accretion onto collapsed objects in the cavity under 



hypotheses similar to our has been recently estimated 
by Hallcr et al. (1996) and turns out to be completely 
negligible. 

Results of model calculations are compared with 
Granat ART-P observations. As discussed in subsec- 
tion 2.1, the diffuse emission detected by Granat in 
the 2.5-8.5 keV band comes from a roughly ellipti- 
cal region of ~ 1.18 deg 2 around Sgr A*. The val- 
ues of the observed total flux reported by Sunyaev, 
Markevitch, & Pavlinsky (1993) are 4.8 ± 1 x 10~ 10 
ergs cm~ 2 s _1 and 4.6 ± 1.3 x 10~ 10 ergs cm~ 2 s _1 
in the two energy bands 2.5-8.5 keV and 8.5-22 keV. 
For the sake of simplicity, we approximate the to- 
tal ellipse with a circle of the same area, centered at 
r = 0. We then calculate the expected X-ray flux in 
the two spectral bands, varying B, v n and Nh- The 
fraction / of neutron stars is always taken equal to 
0.01; obviously fluxes scale linearly with /. Results 
are reported in table 1 . As can be seen, the flux emit- 
ted above 8.5 keV never exceeds a few percent of the 
observed one, so accreting ONSs can not be directly 
responsible for the detected hard X-ray emission. On 
the contrary, the contribution of ONSs in the lower 
energy band can be substantial. For the two values of 
B we have considered, the expected emission ranges 
from a minimum of 0.07 up to ~ 1.6 the Granat flux in 
the 2.5-8.5 keV band. Clearly, this kind of calculation 
is necessarily influenced by the uncertainties in the 
assumed ISM and ONSs distributions. However, al- 
lowing for different values of the ONSs mean velocity 
and of the total column density, computed values are 
always comparable with observational data. In ad- 
dition, while it could be not unreasonable to assume 
higher mean velocities for the NSs population (see e.g. 
Lyne & Lorimer 1994; Haller et al. 1996), the range 
of column densities we have explored corresponds to 
the higher values for the interstellar absorption to- 
ward the Galactic Center reported in the literature 
(see e.g. Thomas et al. 1996). Moreover the fraction 
/ of neutron star has been fixed to a rather conserva- 
tive (low) value. All these considerations suggest that 
the ONSs contribution to the diffuse X-ray emission 
from the Galactic Center can be substantial in the 
low energy component. 

We note that the ONSs integrated luminosity is 
typically ^ 10 38 erg s _1 and that most of the radia- 
tion is in soft X-rays. This may prove of importance 
in the discussion of the heating and dynamics of the 
ISM in the central region of the Galaxy. 

The spectrum of the diffuse emission from both the 
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total 1° x 1.5° ellipse and a smaller circle 30' in di- 
ameter has been obtained by Markcvitch, Sunyaev, 

6 Pavlinsky (1993). For the sake of comparison, we 
have calculated the monochromatic flux from these 
two regions using equation (7). Results are presented 
in figures 1 and 2 for three sets of parameter values. 
Crosses in figures 1 and 2 are the Granat data with 
their error bars. As can be seen, the observed spec- 
tral distribution is flatter than the predicted one, but 
the spectral shape is reasonably reproduced up to 6— 

7 keV. Other physical processes should be invoked to 
explain the harder emission. 

We have compared the radial distribution of the 
emitted flux with available observational data. From 
the central 30' region we expect approximately 34% 
of the total flux. This fraction should be compared 
with the ~ 25% deduced from Granat data. The 
two values are in rough agreement although the radial 
distribution of the diffuse emission predicted by our 
simplified model appears more concentrated toward 
the center than the observed one. Figure 3 shows the 
contour levels of the X-ray flux per square degree, 
normalized to the total, in the 2.5-5 keV band. 

As a final point, we note that accretion onto ONSs 
produces an extremely smooth X-ray source in our 
model. In fact, the neutron star with the highest 
accretion rate contributes only to 3 x 10~ 4 of the total 
X-ray luminosity. 

4. Discussion 

In this paper we have considered the possibility 
that the 2.5-7 keV component of the diffuse X-ray 
source observed in the direction of the Galactic Cen- 
ter is due to the unresolved emission from old neu- 
tron stars accreting the interstellar medium near the 
Galactic Center. Due to the large gas and stellar 
densities, the central region of the Galaxy appears, 
in fact, a very promising site to detect the overall 
emission from accreting ONSs. We have shown that 
a sizable fraction, possibly most, of the 2.5-7 keV 
emission observed in the Galactic Center may be ex- 
plained in terms of an old neutron star population 
accreting the dense ISM. Because of the intrinsic un- 
certainties in our assumptions, mainly in the value 
of the neutron star relic magnetic field and in the 
star velocity distribution, it is not possible to assess 
firmly that ONSs are responsible for the extended X- 
ray emission in this band. Moreover, accretion itself 
onto magnetized neutron stars could be questionable. 



Although the propeller effect is not a seriuos threat 
if B ~ 10 9 G (see Blaes, & Madau 1993; Treves, 
Colpi, & Lipunov 1993), preheating of the infalling 
gas due to the emitted radiation may inhibit accre- 
tion, as recently pointed out by Blaes, Warren, & 
Madau (1995). Modulo these caveats, our calcula- 
tions show that the observed intensity, spectral shape 
and flux radial dependence are substantially well re- 
produced for different values of the free parameters of 
the model. Above ~ 6 — 7 keV the integrated ONSs 
spectrum gives only a marginal contribution to the 
observed flux. 

The calculated continuum is rather flat and then 
drops sharply above ~ 6 keV, so we do not expect 
that ONSs emission could provide an important pho- 
toionization source to produce ultimately iron lines. 
However, our synthetic spectra were computed con- 
sidering only pure hydrogen atmospheres around ac- 
creting neutron stars. A definite assessment about 
the presence of Ka lines in the emerging spectrum 
would require the extension of our model to include 
line processes and a more realistic chemical composi- 
tion for the accreting gas. If accretion onto neutron 
stars with a magnetic field as high as ~ 10 12 G is 
possible (which is dubious, see, e.g., Treves, Colpi, 
& Lipunov 1993), about 10 % of the flux is emitted 
in the cyclotron line at Eb = II.6B12 keV (Nelson et 
al. 1995) and the excitation of the Fe lines could be 
substantial. 

Sunyaev, Markevitch, & Pavlinsky (1993) noticed 
some correlation between the distribution of giant 
molecular clouds and the X-ray brightness and inter- 
preted it as a signature of the reprocessing of medium 
energy photons by the clouds themselves. Such a cor- 
relation is naturally explained in our scenario, since 
GMCs provide regions of very high density in which 
the accretion rate can be substantially higher than 
our present estimate, based on n = 100 cm -3 . This 
implies also harder spectra, so that the emission of 
hard photons, ~ 10 keV, is possible for ONSs accret- 
ing in clouds with n ~ 10 4 cm~ 3 . On this regard, see 
the discussion of Campana & Mercghctti (1993) on 
the source 1740.7-2942. 

The correlation between X-ray emission and the 
ISM distribution may become the basic probe to test 
the real relevance of ONSs accretion in explaining the 
diffuse X-ray source in the Galactic Center. AXAF 
and JET-X will have the required sensitivity and 
space resolution in medium energy X-rays to shed 
light on this important issue. 
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Table 1 

Computed fluxes for different parameter values 





N H 


S 


^[2.5-8.5] 


^[8.5-22] 


km s _1 


10 23 cm- 2 


G 


10~ 10 erg cm -2 s -1 


10~ 10 erg cm" 


75 


2 


10 9 


7.11 


0.071 


75 


3 


10 9 


3.90 


0.064 


75 


4 


10 9 


2.44 


0.057 


100 


2 


10 9 


3.18 


0.031 


100 


2.5 


10 9 


2.30 


0.030 


100 


3 


10 9 


1.74 


0.028 


100 


3.5 


10 9 


1.36 


0.027 


100 


4 


10 9 


1.09 


0.025 


125 


2 


10 9 


1.68 


0.017 


125 


3 


10 9 


0.92 


0.015 


125 


4 


10 9 


0.57 


0.013 


150 


2 


10 9 


0.99 


0.010 


100 


2 


10 10 


7.80 


0.151 


100 


3 


10 10 


4.61 


0.135 


100 


4 


10 10 


3.03 


0.121 


200 


2 


10 10 


1.05 


0.020 


300 


2 


10 10 


0.32 


0.006 



Note.— / = 10 2 , a = 1 pc, z max = 300 pc 
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Fig. 1— Granat ART-P data (Markevitch, Sun- 
yaev, & Pavlinski 1995, crosses) and calculated X-ray 
spectra from a region of 1.18 deg 2 for: a) / = 0.01, 
v = 75 km/s and Nh = 3 x 10 23 cm~ 2 (solid line); 
b) / = 0.02, v = 100 km/s and N H = 3 x 10 23 
cm -2 (dashed line); c) / = 0.015, v = 100 km/s and 
N H = 2.5 x 10 23 cm" 2 (dash-dotted line). For all 
models B = 10 9 G. 

Fig. 2. — Same as in figure 1 for the central region 
30' in diameter. 

Fig. 3. — Contour levels of the X-ray flux/deg 2 in the 
2.5-5 keV band. Parameters values are B = 10 9 G, 
/ = 0.01, v a = 75 km/s and N H = 3 x 10 23 cur 2 . 
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